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(N ; Abstract 

The formalism of Supersymmetric Quantum Mechanics supplies a trial wave function 
to be used in the Variational Method. The screened Coulomb potential is analysed 
within this approach. Numerical and exact results for energy eigenvalues are com- 



O ! pared. 
(N 

q ■ I. Introduction 

O ' 

The screened Coulomb potential has been used in several branches of Physics, for 
instance, in nuclear physics (as the name of Yukawa potential), in plasma and in the 
q \ study of electrolytic solution properties (Debye-Huckel potential). The Schrodinger 
equation for this potential is not exactly solvable and exact numerical, [p]], || and 
approximative, ||-[f| methods have been applied to obtain the energy eigenvalues, 
including variational calculations, ||-||. 

More recently, a new methodology based within the variational method associ- 



ated to supersymmetric quantum mechanics formalism has been introduced, Pl-fllT 



References and |T(| introduce a scheme based in the hierarchy of Hamiltonians; 



it permits the evaluation of excited states for one-dimensional systems. In reference 
[ll| an ansatz for the superpotential which is related to the trial wave function is 
proposed. The new approach showed to be useful to get answers when applied to 
atomic systems, Jll]|-[0. 

In this letter, the variational energy eigenvalues for the static screened Coulomb 
potential in three dimensions are determined using the variational method using a 
trial wave function induced by Supersymmetric Quantum Mechanics, SQM. 

In the approach followed here the first step taken is to look for an effective potential 
similar to the original screened Coulomb potential. Inspired by SQM, an ansatz is 
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made to the superpotential which determines the variational (trial) wavefunction 
through the algebraic approach of SQM. 

Our system is three dimensional and in this case it is possible to determine the 
variational eigenfunctions for each value of angular momentum I. The first eigenfunc- 
tion, obtained from direct factorization of the effective Hamiltonian, corresponds to 
the minimum energy for each /. 

This new methodology has been successfully applied to other atomic systems such 
as the Hulthen, |Q), and the Morse, [12|], potentials. Here it is applied to the screened 



Coulomb potential. 

In what follows, we briefly introduce, for completeness, the SQM scheme, then 
introduce the variational method and show our results. 

II. Supersymmetric Quantum Mechanics 



In SQM, fIQ|-[n|, for iV = 2 we have two nilpotent operators, Q and Q + , that 
satisfying the algebra 

{Q,Q + } = H SS ;Q 2 = Q+ 2 = 0, (1) 
where Hss is the supersymmetric Hamiltonian. This algebra can be realized as 

where are bosonic operators. With this realization the supersymmetric Hamilto- 
nian Hss is then given by 

Hss = ( A 4 A-A+ ) = ( H H- ) ' (3) 

where if ± are called supersymmetric partner Hamiltonians and share the same spec- 
tra, apart from the nondegenerate ground state. Using the super-algebra, a given 
Hamiltonian can be factorized in terms of the bosonic operators. In H = c = 1 units, 
it is given by 

Hi = -\§r 2 +V l {r) = AtA- l +E$ ) (4) 

where is the lowest eigenvalue and the function Vi(r) includes the barrier poten- 
tial term. The bosonic operators are defined by 

where the superpotential W±(r) satisfies the Riccati equation 

Wl - W[ = 2Vx(r) - 2E^ ] (6) 
which is a consequence of the factorization of the Hamiltonian H\. 
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The eigenfunction for the lowest state is related to the superpotential W\ by 



*W( r ) = Nexp(- f Wx{f)df). (7) 
J o 

Now it is possible to construct the supersymmetric partner Hamiltonian 



H 2 = A', At + £f = -~ + \{Wl + W[) + 4 1} 



If one factorizes H 2 in terms of a new pair of bosonic operators, A 2 one gets, 

H 2 = AiA~ 2 + £f = + i(W 2 2 - W*) + (9) 

where Eq is the lowest eigenvalue of H 2 and W2 satisfy the Riccati equation, 

Wl -W 2 = 2V 2 (r) - 2E ( 2) . (10) 

Thus a whole hierarchy of Hamiltonians can be constructed, with simple relations 
connecting the eigenvalues and eigenf unctions of the n-members, |fl3| , |15[ . 



Thus, the formalism of SQM allows us to evaluate the ground state eigenfunc- 
tion from the knowledge of the superpotential W(r), satisfying the Riccati equation, 
eq.(||). However, since the potential is not exactly solvable, the purpose is to propose 
an ansatz for the superpotential and, based in the superalgebra information, we eval- 
uate a trial wavefunction that minimizes the expectation value of the energy. The 
energy eigenvalues pursued are evaluated by minimizing the energy expectation value 
with respect to a free parameter introduced by the ansatz. 

III. Trial wavefunction for the variational calculation 

The screened Coulomb potential is given, in atomic units, by 

p-Sr 

V SC = — (11) 

r 

where 6 is the screened length. The associated radial Schrodinger equation includes 
the potential barrier term and it is given by 

/ 1 d 2 e~ Sr l(l + l)\ T ^ T 

h v ' \I> = E^f (12) 

\ 2dr 2 r 2r 2 j K } 

where the unit length is h 2 /me 2 and the energy unit is eo = —me^/h 2 . 

In order to determine an effective potential similar to the potential in the Hamil- 
tonian (|T2"|), that is the screened Coulomb potential plus the potential barrier term, 
the following ansatz to the superpotential is suggested 

5e~ Sr 1 5 

W(r) = -(l + l)— —— + ——---. (13) 
v ; v ' 1 - e~ 5r Z + 1 2 v ; 
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Substituting it into (0), one gets 



¥ ( r ) = (1 - e- Sr ) l+1 e' { oTT)'^ r . (14) 

Assuming that the radial trial wave function is given by fll4]) , replacing S by the 
variational parameter /i, i.e., 

* M (r) = (1 - e -^y+i e -(izri7-f> (15) 
the variational energy is given by 

f _ jTMlhi£^izi±MlM^ run 

Thus, minimizing this energy with respect to the variational parameter fj, one obtains 
the best estimate for the energy of the screened Coulomb potential. 

As our potential is not exactly solvable, the superpotential given by eq.(O) does 
not satisfy the Riccati equation (|6|) but it does satisfy it for an effective potential 
instead, V e ff 

W? - Wi 

V ef f(r) = 1 2 1 + E(fl) (17) 

where W\ = W±(S = fi) is given by eq. (|l3D and is the parameter that minimises the 
energy expectation value, (|T^|). It is given by 

Veffir) = + 2 -7T3^ + 2 ( 7TT - 2 } + (18) 



e 



where 5 = jl that minimises energy expectation value. One observes that for small 
values of 5 the first term is similar to the potential (|TTD and the last is approximately 
the potential barrier term. This observation allows us to conclude that the superpo- 
tential ( |13D can be used to analyse the three dimensional screened Coulomb potential 
variationally through the trial wavefunction (|I4|). 



IV. Results 



For I = the effective potential (|j|) becomes identical to the Hulthen potential. 
Thus, the results presented in Table 1 coincide with those of ref.0, where the Hulthen 
potential eigenfunctions are directly used in the variational calculation. The deviation 
on the fifth decimal algarism can be attributed to the accuracy of the numerical 
calculation. 
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State Is 


SfTPPTliriP' A 

kJ V_yJ_ VjV^J.J.±±J. fcl \J 


SOM Variational 

U V/1V1 V (Ail. ICAi U J.V_y J. J.CAX 


Variational (T~?pf 7) 


T^/va rt IMnmPTiral 

J /vVOjVj U J. 1 1. 1. 1 J. J. V - i_ 1 V. < A.l 


0.001 


-0 49902 


-0 49900 




0.002 


49802 


-0 49800 

vy • * ./ vj vy vy 


-0.4980 


0.005 


-0 49504 

vy • *../ vy vy T: 


-0 49502 

vy • * ./ vy vy 


-0.4950 


0.010 


-0 49009 


-0 49007 


-0.4901 


0.02 


-0 48031 

Vy . j: VJ Vy VJ X 


-0 48030 

W . j: VJ Vy VJ Vy 


-0.4803 


0.025 


-047548 


-0.47546 


-0.4755 


0.03 


-0 47068 


-0 47066 

vy • i vy vy vy 




0.04 


-0.46119 


-0.46117 




0.05 


-0 45180 

V 7 . 1UU 


-0 45182 

Vy • t^y ilju 


-0.4518 


0.06 


-0 44259 


-0 44260 




0.07 


-0 43351 


-0 43352 




0.08 


-0 42456 

Vy • JI ^.y V./ 


-0.42457 




0.09 


-0.41574 


-0.41575 




0.10 


-0 40705 


-0 40706 

Vy . i Vy 1 Vy V ' 


-0.4071 


0.20 


-0 32681 


-0 32681 


-0 3268 


0.25 


-0 29092 


-0 29092 


-0 2909 


0.30 


-0 25764 


-0 25763 




0.40 


-0 19842 


-0 19836 

Vy • J-(7UUU 




0.50 


-0 14806 

Vy . ±TUUU 


-0 14808 

Vy . X^vyVyVJ 


-0.1481 


0.60 


-0106077 


-0 10608 

vy • -L vy vy vy vy 




0.70 


-0.07175 


-0.07174 




0.80 


-0.04459 


-0.04459 




0.90 


-0.02420 


-0.02418 




1.00 


-0.01026 


-0.01016 


-0.01029 


1.05 


-0.00568 


-0.00544 





Table 1. Energy eigenvalues as function of the screening parameters 5 for Is 
state (/ = 0), in rydberg units of energy. Comparison is make with variational and 
exact numerical results from [|IJ and M. 

The results become more interesting for 1^0. In this case the effective potential 
differs from the Hulthen potential. Table 2 shows the results for 2p (1 = 1), 3d (I = 2) 
and 4/ (I = 3) energy levels. Also given in this table are the correponding numerical 
results PJ. 
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2] 


3 


3d 


4f 


5 


Variational 


Numerical 

_L 1 I L111V.. 1 X V < A 1. 


"Variational 

V CX J. ± Oj U 1 V_y A J- C(j 1 


TMnmpripal 

X 1 I. i. -L -L -L V yi_ 1 V < A 1 


Va "national 


IMirmprira 1 


0.001 


-0.2480 


-0.2480 


-0 10910 


-0 10910 


-0 06051 


-0 06052 


0.005 










-0.52930 


-0.05294 


0.010 


-0.2305 


-0.2305 


-0.09212 


-0.09212 


-0.04419 


-0.04420 


0.020 


-0.2119 


-0.2119 


-0.07503 


-0.07503 


-0.02897 


-0.02898 


0.025 


-0.2030 


-0.2030 


-0.06714 


-0.06715 






0.050 


-0.1615 


-0.1615 


-0.03374 


-0.03383 






0.100 


-0.09289 


-0.09307 











Table 2. Energy eigenvalues as function of the screening parameters 5 for 2p 
(/ = 1), 3d (I = 2) and 4/ (I = 3) states, in rydberg units of energy. Variational 
values obtained by the trial function (|i~5"D are compared with exact numerical results 
obtained from reference |1[, (see also |J). 

V. Conclusions 

We have proposed trial wavefunctions to be used in the variational calculation in 
order to determine the energy eigenvalues of the screened Coulomb potential. These 
functions were induced from supersymmetric quantum mechanics formalism. The 
approach consists of making an ansatz in the superpotential which satifies the Riccati 
equation by an effective potential. The trial wavefunctions were then determined 
from this superpotential through the superalgebra. 

For / = the effective potential obtained is identical to the Hulthen potential. 
However for / ^ the potential has a new structure. The trial wavefunctions sug- 
gested for this case are different from those proposed in references [§fl-|§. Within 
our framework the energy eigenvalue for each value of I is obtained using the same 



function ([14]). 

In terms of the hierarchy of Hamiltonians, we obtained the first member for each 
value of /. Other members can be determined from the usual approach in supersym- 



metric quantum mechanics, [ l5| . 



One observes that the results obtained are in very good agreement to those found 
in the literature. The results are better for small values of the parameter 5. This 
observation is justified by the fact that for small values of S the effective potential is 
more similar to the original potential than for higher values of 5. 

We stress that even though the problem has been attacked by different methods 
our new methodology is very simple to supply accurate results. We believe that other 
applications to atomic physics problems can be made by this new method. 
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